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Abstract 

The energy spectra of gravitational waves (GW) produced in quintessential inflationary 
models increase in frequency and exhibit a sharp spike around 170 GHz where the associated 
^ ■ fraction of critical energy density today stored in relic gravitons is of the order of 10 -6 . We 
contrast our findings with the spectra of ordinary inflationary models and we comment about 
Q\ ■ possible detetction strategies of the spike. 
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In ordinary inflationary models, Oqw (the present fraction of critical energy density 
stored in relic gravitons) is notoriously quite small. In fact Oqw is either flat of decreasing as a 
function of the present frequency. Therefore, the COBE bound, applied at the frequency scale 
of the present horizon (i.e. vq ~ 1.1 x 1CT 18 h Hz, 0.5 < h < 1), demands h% Qqw < 6.9 x 
10~ n [[[]. Since the energy spectrum decreases sharply as v~ 2 between v and the decoupling 
frequency, we can further argue that for u > 10~ 16 Hz, /iq^qw cannot exceed 10~ 14 . This 
conclusion can be evaded provided the inflationary phase is not followed immediately by the 
radiation dominated epoch but rather by an expanding phase driven by an effective source 
whose equation of state is stiffer than radiation. Indeed a stochastic background of relic 
gravitons can be produced, with different spectra, in any variation of the expansion rate of 
the Universe @. 

Recently Peebles and Vilenkin discussed a model where the occurrence of a stiff (post- 
inflationary) phase can be dynamically realized ||. As previously argued [|J the graviton 
energy spectra in this class of models must be increasing as a function of the present fre- 
quency. One of the motivations of || stems from a recent set of observations which seem 
to imply that Q m (the present density parameter in baryonic plus dark matter) should be 
significantly smaller than one and probably of the order of 0.3. If the Universe is flat, the 
relation between luminosity and red-shift observed for Type la supernovae || hints that 
the missing energy might be stored in a fluid with negative pressure acting as an effective 
(time dependent) cosmological term whose magnitude should be of the order of 10~ 47 GeV 4 , 
too small if compared with the cosmological constant arising from electroweak spontaneous 
symmetry bereaking (which would contribute with (250 GeV) 4 ). A complementary way of 
thinking is that the missing energy could come from a dynamical scalar (the quintessence 
H field) whose potential is unbounded from below [0. The starting point of || is that 
could be identified with the inflaton and, as a consequence of this identification the effective 
potential of will have to inflate for < and it will be unbounded from below for > 0. 
As an example we could take 

AAf 8 

V(0) = A(0 4 + M 4 ), for0<O, and V{<f>) = , 4 + M 4 > for > 0. (1) 

where, if we want the present energy density in to be comparable with (but less then) 
the total (present) energy density we have to require M ~ 10 6 GeV. Any other inflationary 
potential can be used for < 0. 

For a long period after the end of inflation the kinetic term of will dominate the stress 
tensor and, therefore, the effective fluid driving the geometry will have a speed of sound equal 
to the speed of light. The energy spectra of the relic gravitons will then be blue, namely they 
increase with frequency with a power wich depends, in general, upon the precise equation 
of state Not only gravitons are parametrically amplified in this class of models but 
also any other (non conformally coupled) scalar degree of freedom |8|, |9fl. During the stiff 
phase the energy density of the produced fluctuations red-shifts more slowly than the energy 
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density of the background, and, at some moment the energy density of the produced quanta 
will become dominant triggering the reheating of the Universe |§ . GW and inflaton quanta 
(equaivalent to 3 degrees of freedom) are unable to reheat the Universe on their own: their 
spectra, are non thermal j|] and cannot thermalize below the Planck scale. If iV s minimally 
coupled scalar field are present they can reheat the Universe with a thermal distribution 
since their energy spectra, amplified because of the transition from the inflationary to the 
stiff phase, can thermalize thanks to non-graviational (i.e. gauge) interactions which get 
to local thermal equilibrium well below the Planck energy scale. The Universe will become 
eventually dominated by radiation. This will occur at a temperature which is a function of 
Hi, the curvature scale at the end of inflation, and of N s : 

T r = (^)R 3/4m p - !0 3 ^ s 3/4 GeV, R = N s R i: R, ~ 1(T 2 . (2) 

If we do not fine-tune Hi to be much smaller than 10~ 7 in Planck units, T r is typically a bit 
larger than 1 TeV. Ri is the fractional contribution of each (minimally coupled) scalar degrees 
of freedom to the energy density of the produced quanta right after the end of inflation. 

In this letter we are interested in the calculation of the energy spectra of the pure (trans- 
verse and traceless) tensor modes of the geometry 

g^{x,r)) = a 2 (r])[r} f t U + hi ll/ (x,r))), with h fl0 = 0, = °> h^ = 0, (3) 

where rj^ is the usual Minkovski metric and V M is the covariant derivative associated with 
the (conformally flat) background geometry. We will focus our attention on the hard branch 
of the spectrum namely on those tensor modes which left the horizon before the end of 
inflation and re-entered during the stiff phase. Since GW only couple to the curvature and 
not to the matter sources (which can only support scalar inhomogeneities) the spectrum will 
be fully determined by the scale factors whose evolution reads, in conformal time, 
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for r\ < —rji, and, a s (rj) = J - ~*~ — — , for — 771 < 77 < r\ r (4) 



where 771 = (aiHi)' 1 and H r = (a^) -1 is the curvature scale at the temperature T r when 
the radiation phase commences. Notice that in Eq. (f|) the scale factors and their first 
derivatives (with respect to the conformal time if) are continuous in —771. 

The mode function associated with the two polarization of stochastically distributed GW 
obeys the (Schroedinger-like) equation 



k 2 - = 0, tfj = ah, ' = 4- (5) 

or] 



which has to be solved in each of the two temporal regions defined by Eq. (|j). Given the form 
of a" /a in the case of Eq. (HI), ip will be a linear combination of Bessel functions, oscillating 
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for k 2 ^> \a"/a\ but parametrically amplified in the opposite limit (i.e. k < \a"/a\): 
if) l {k^) = ^=^ 2 \x), p =] jle-^ +1 \ ri<^ 

MKv) = ^§f[s*A + (k)H! ) 2 \y) + sA4k)Hi ) 1 \y)}, s = J^j*, - Vl < V < Vr , (6) 
V 2k V l 



where x = kr\ and y = k(r) + p and s guarantee that the large argument limit of the 
Hankel functions Hj) 1 ' 2 ^ is exactly the one required by the quantum mechanical normalization 
(namely e ±lkri / y/k) . In the case of a pure de Sitter phase v = 1.5 but corrections (of few 
percents) can arise if the slow-rolling corrections are taken into account [ |T0|j . 

The graviton energy density per logarithmic interval of frequency will then be given by 

pu> = d f, GW = n(u) = \A_(u)\ 2 , oj = - = 2irv, (7) 

where u is the physical wavenumber and v the physical frequency. Because of the continuity 
of a(rj) and a'(rj), the two mixing coefficients A±(k) can be fixed by the two conditions 
obtained matching if) and if)' in rj = —rji with the result that 

Q 

A(k) ~ — 2 v -h-i n{2u+1) r(u)x^ u \nx 1 , valid for x x < 1. (8) 

Notice that for X\ > 1 the mixing of the modes is exponentially suppressed and the ultra- 
violet divergence is avoided ||, ffl. Inserting Eq. (H) into Eq. (H) we get the hard branch of 
the relic graviton energy spectrum (in critical units) 

which is defined, at the present time 770, between the two frequencies 

u J Vo ) = 3.58 BJ ( -\r) f^V /3 mHz, and = 358 iT*(— GHz, (10) 

MO 14 / V # th / V g t h 



^■G\v(^,Vo) = — = fi 7 (*o) £ (jf~) (~~) lll2 (^)' W r <W<Wl, (9) 



^dec\ 1/3 n f . v _ P 7 (* ) _ 00^ T 4 - ., 



where 

V 0m / Pc(*o) 30 if z M|, 

fi 7 (to) is the fraction of critical energy density in the form of radiation at the present obser- 
vation time; go = 2, T = 2.73 K; g^cc = 3.36 and g t h = 106.75 are, respectively, the number 
of (massless) spin degrees of freedom at decoupling and at thermalization. The dependence 
upon the number of relativistic degrees of freedom in f^cw occurs since, unlike gravitons, 
matter thermalizes and then the ratio between the pew an d p c is only approximately constant 
in the radiation dominated phase. 

By taking H\/Mp = vA < 10~ 7 the spectrum satisfies the COBE bound [|IJ. Since the 
spectral energy density increases sharply in the hard branch the most relevant constraints 
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will not come from large scales (as in the case of ordinary inflationary models) but from 
short distance physics and, in particular, from big-bang nucleosynthesis (BBN). In order to 
prevent the Universe from expanding too fast at BBN we have to demand 
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d\nun GW (u,t n ) < — (N u -3) 



Pc(*n) 



;i2) 



where t n is the nucleosynthesis time. Since the number of massless neutrinos N u cannot 
exceed, in the homogeneous and isotropic BBN scenario is bounded, 3.4 we have that the 
nucleosynthesis bound implies 



fjn 



1/3 



— — ) <0.07, 



(13) 



where the factor of 3 counts the two polarizations of the gravitons but also the quanta 
associated with the inflaton and g n = 10.75 is the number of spin degrees of freedom at 
t n . From Eq. fll~3|), N s > 19.9 as it can occur, for instance, in the minimal supersymmetric 
standard model (MSSM) where N s = 104 but not in the Minimal Standard model where 
there is only one Higgs doublet with two (complex) scalars and N s = 4 ]3[]. 

From Eqs.fl9D-(|T3"D, assuming, for instance, N s = 21 the present coordinates of the relic 
graviton spike, providing an overall normalization of the whole spectrum, will be 



1/1(770) = 170 GHz, /igflGw(fi, Vo) = 0-8 x 10" 



(14) 



eight order of magnitude larger than the signal provided by ordinary inflationary models 
Ll, 12]. An increase in N s decreases the height of the spike. The decrease is quite mild 
since, from Eq. @ we can deduce that, at the spike, f^u^r/ci) oc iV~ 3 / 4 . An increase in 
iV s makes narrower the peak structure associated with the spike. In fact z/ r (^ ) oc N^ 4 gets 
larger for larger N s whereas ^i(?7o) oc A"" 1 / 4 gets pushed towards more infra-red values of 
the spectrum. 

In Fig. [I] we report plot the spike computed from Eq. (H) by taking into account the 
bound of Eq. (|T3|). A decrease in the curvature scale at the end of inflation does not affect the 
spike and /or the maximal frequency of the spectrum since v\ does not depend on Hi/Mp. In 
the case of ordinary inflationary models z^i (?7o) = 100^J Hi/Mp GHz and, for Hi < 10 _7 Mp, 
Pi (rjo) can be, at most, 0.1 GHz. In our case, by decreasing Hi, ^i(?7o) does not change but 
Vr{vo) decreases making the peak broader (in Fig. [I] a decrease in v r (vo) shifts the starting 
point of the hard branch to the left keeping fixed the position of the spike). For instance 
if Hi = 10~ 7 Mp (as assumed in the plot of Fig. |1|) the spike is localized according to Eq. 
(|T4]) and z/ r (?]o) = 170 mHz giving a range of fourteen orders of magnitude where the energy 
density increases as u\nu. 

In the next few years various interferometric detectors like LIGO |]13" |, VIRGO ||14|| , GEO- 



600 |15j will come in operation. The spectral densities of the the noise are peculiar of each 
detector but they are all defined between 1 Hz and 1 kHz with a maximal sensitivity around 
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Figure 1: We illustrate the energy spectrum of the relic GW from quintessential inflation as 
a function of the physical wavenumber. We normalized the height of the spike appearing in 
the hard branch (see Eq. @) to be compatible with the BBN bound. With black spots we 
indicated the COBE and millisecond pulsar timing bound and the approximate Virgo/LIGO 
and advanced LIGO sensitivity. The shaded area does corespond to the region where the 
spike is above the signal provided by ordinary inflationary models. For completeness we 
indicated also the soft and semi-hard branches of the spectrum whose detailed calculation 
can be found in ||10 



0.1 kHz. In this frequency range the spectral density of the signal, Shiy) can be related to 
the energy density: 

47T 2 

M G w(v,Vo) = ^772V 3 S h (v,T] ). (15) 

Using now Eq. (||) into Eq. (|15|) we have 



S h (u;, m )=CR-i (^)- 1 ^ 7 (t ) (-Y 2 ln 2 f-) Hz" 1 , (16) 

with C = 6.5 x 10" 73 hi For u ~ 0.1 kHz, S h ~ lO^-lO" 53 sec. For u ~ 0.01 kHz, 
Sh ~ 10~ 50 -10~ 51 sec. The spectral density of our signal should be carefully compared 
with the spectral density of the noise. Our signal is too weak to be interesting for the first 
generation of interferometers. The sensitivity which is closer to the signal of quintessential 
inflationary models correposnds to the case of the two upgraded LIGO detectors. 

Let us estimate the strength of our background for a frequency of the order of 0.1 kHz 
-1 kHz. Let us assume that the energy density of the stochastic background is the maximal 
compatible with the nucleosynthesis indications. As a function of N s , the GW energy density 
(in critical units) at a frequency uj ~ 0.1-1 kHz is then 

JW(i/j, 770) h 2 = 2.29 10~ 15 iV7 3/4 [-19.7 + 0.25 In N s ] 2 , vi = 0.1 kHz, (17) 
n GW (i/ ; , % ) hi = 2.29 10~ 14 iV7 3/4 [-17.4 + 0.25 In iV s ] 2 , Vl = 1 kHz. (18) 
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Suppose then that we correlate the two LIGO detectors for a period r = 3 months. Then, 



the signal to noise ratio (squared) can be expressed as |L6 

S\ 2 9# 4 _ r...7>)°GwMo) 



N) SOvWo d \*sV(v)S%\vy (19) 

where is the overlap function accounting for the difference in location and orientation 
of the two detectors. For detectors very close and parallel, 7(f) = 1. For the two LIGO 
detectors [] and for other pairs of detectors 7(1/) can be computed flIE| . Sfy 2 \v) are the noise 
spectral densities of LIGO-WA and LIGO-LA and since the two detectors are supposed to 
be identical we will have that S$(u) = S$(u). In order to detect a stochastic background 
with 90% confidence we have to demand S/N>,1.65. For an estimate of S/N we need to 
evaluate numerically the integral appearing in Eq. (|19|) where the theoretical information 
comes from HgWj given, in our case, by Eq. (|9|) The experimental information is encoded in 
the noise spectral densities of the LIGO detectors which are not of public availability. We are 
not aware of any calculation of the sensitivity of the LIGO detectors for an energy spectrum 
whose frequency behavior is the one of Eq. In the case of a flat energy spectrum the 
S/N has been computed [I7[] and we have that the minimum Oqw detectable in r = 4 



months is given, with 90 % confidence, by ficw (i/j,7fo) = 5 x I0~ 6 h 2 (for the initial LIGO 
detectors) and by Qqw^vj, rj ) = 5 x 10 -11 /iQ 2 (for the advanced LIGO detectors) JT71. In 
Fig. we compared our signal given, at the interferometers frequency vj y by Eqs. ( pf ) and 
(|18|) with the sensitivity of the advannced LIGO project to a flat spectrum. For the allowed 
range of variation of N s our signal lies always below (of roughly 1.5 orders of magnitude) 
the predicted sensitivity for the detection, by the advanced LIGO, of an energy density with 
flat slope. The main uncertainty in this analysis is however the spectral behavior of the 
sensitivity for a spectrum which, unlike the one used for comparison, is not flat. It might be 
quite interesting to perform accuarately the calculation of the S/N in order to see which is 
the precise sensitivity of the LIGO detectors to a spectral energy density as large as 10~ 12 
and rising as {y jv r } In [y/v\) in a frequency range 1 Hz-1 kHz. 

On top of the interferometric and resonant bar detectors electromagnetic detectors and, in 
particular microwave cavities could be employed, in the future, in order to detect background 
of relic gravitons coming from quintessential inflation. In fact the the nucleosynthesis bound 
is almost saturated for frequencies of the order of v\ = 358 x i? -1 / 4 GHz. Microwave 
cavities can be used as gravitational waves detectors in the GHz frequency range There 
were published results reporting the construction of such a detector fL9| . In this prototype 
i/Gw = 10 GHz and the sensitivity to fractional deformations 5x/x was the order of 10~ 17 
using an integration time At ~ 10 3 sec. There, are at the moment, no operating prototypes 
of these detectors and so it is difficult to evaluate their sensitivity. The example we quoted 
refers to 1978. We think that possible improvements especially in the quality factors of 



2 One LIGO detector (LIGO-WA) is being built in Hanford (near Washington) the other detector (LIGO- 
LA) is under construction in Livingston (Lousiana). 
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Figure 2: We illustrate the signal of the quitessential graviton background at the frequency 
of the interferometers. The two thin lines correspond to the signal at uj = 0.1 kHz and 
vj = 1 kHz as a function of N s according to Eqs. (17) and (18). The region between the 
thin lines represent approximately our signal and the full thick line represents the sensitivity 
of the advanced LIGO detectors to a stochastic background with flat energy spectrum. In 
order to be detected our signal should lie in the dashed area. We see that for the allowed 
range of variation of N s the signal is always smaller than the sensitivity. This comparison 
is only illustrative and not completely correct. In fact, the sensitivity our specific energy 
spectrum (increasing as ulnu), is not expected to be exactly equal to the sensitivity to a flat 
flaw- The precise sensitivity is, in principle, computable and it requires, according to Eq. 
(19), the knowledge of the spectral density of the noises which are not publically available. 
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the resonators can be envisaged. In spite of the fact that improvements can be foreseen we 
can notice immediately that, perhaps, to look in the highest possible frequency range of our 
model is not the best thing to do. In fact we can argue that in order to detect a signal of 
the order of /iq^gw ~ 10~ 6 at a frequency of 1 GHz, we would need a spectral density of 
the noise smaller than the one of the signal, which, from Eq. (|T5| ) turns out to be 

/kHz\ 3 

S h {u,Vo)^ x 1(T 52 (--) sec, (20) 
corresponding to a sensitivity to fractional deformations of the order of 10 -30 . Moreover, 



as stressed in [20 and already noticed in [19] the thermal noise is one of the fundamental 



source of limitation of the sensitivity of these detectors. An interesting strategy could be to 
decrease the operating frequency range of the device by going at frequencies of the order of 
1 MHz. 

The common lore is that inflationary models cannot give rise to large energy densities 
stored in relic gravitons. We showed that this conclusion is in fact evaded if, as in the 
case of quintessential inflation, a stiff phase follows the inflating epoch. The resulting signal 
can then be eight orders of magnitude larger than the one obtained in the case of a direct 
transition from an inflating epoch to a radiation phase. At the LIGO frequency our signal 
is just below the advanced LIGO sensitivity to flat spectra. Concerning the detectability of 
our background two final comments are in order. The LIGO-LA/LIGO-WA sensitivity to 
our specific spectra has not been computed and we wonder if this could be perhaps done 
in the future. The GHz region (where our signal is maximal) should be carefully explored 
perhaps with the use of electromagnetic detectors. 

I would like to express my gratitude to A. Vilenkin for very useful comments and sug- 
gestions which stimulated the present investigation. 
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